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ABSTRACT 

This paper presents the latest results from experiment E787, at Brookhaven 
National Laboratory, on K + -^>-k + vv and radiative K + decays. The result for 
K + -^-k + vv uses data collected in runs taken during 1995, 1996 and 1997. In 
addition, we discuss plans for future measurements of K + -^tt + i/V . 



1 Theoretical Background 

1.1 K+^TT+VV 

This decay mode is a very clean probe of the dynamics of the Cabbibo-Kobayashi- 
Maskawa (CKM) mixing matrix. It is a AS = 1 process which proceeds via the 
box diagram and electroweak penguin diagrams 0) . Since the top quark is so 
much more massive than the other up-type quarks, its contribution dominates 



the decay rate. In addition, long distance effects are small and the hadronic 
matrix element can be obtained from the semi-leptonic decay K + — > Tt a e + v e 
(K e 3). As a result, the theoretical uncertainty in understanding the decay rate 
is rather low (theoretical error is « 7%, most of it is due to the charm quark 
contribution). The rate can be expressed (in the Standard Model) as, 

B(K+^»v) = l +a \ B } K ^ - ^\X t X t + X%\" (1) 
2tt 2 sin 4 9w\VusV 

where, k + is the isospin correction, A t = V ts V t * d , X c = V cs V* d , and X t ^ are 
Inami-Lim functions 0) , and the sum is over the three neutrino flavours. 
Using measurements from the K and B systems, to determine parameters of 
the CKM matrix, which are used as inputs, one can predict, 

0.5 x 1(T 10 < B(K+^w + i>V ) <1.2xl(r 10 

This range is due to the uncertainties in the input parameters. Since the 

top quark dominates the proceedings, this decay mode provides a very clean 

measurement of Vtd- This measurement is complementary to the measurement 
n — o 

of V t d from B d B d mixing. A difference in the value of Vtd extracted from these 
two sources, could suggest new physics, since any new phenomena would, in 
general, affect the K and B systems differently. 

1.2 Radiative K + decays 

Radiative K + decays come from two sources, (a) the "pedestrian" Inner 
Bremsstrahlung decay, where a charged decay daughter emits a photon, and 
(b) the more interesting structure dependent (SD) part. In this paper, we are 
concerned with the SD contribution. 

The latter source is a good testing ground for Chiral Perturbation theory 
(xPT), and is also important for measuring long distance contributions to other 
decays of interest, e.g., K Q L — » is a background to K\\ — ► ir°l + l~. 

For K + ^> n + v^ , form factors for the SD part, Fa and Fy are predicted 
by xPT. K + ^7r + 7r°7 is interesting in its own right; decay rates for K + , Kl 
and Ks could be similar, even though the decay rates for the non-photonic 
final state (mr) are very dissimilar, also, some SUSY models predict that K+ 
and K~ decay rates are unequal, leading to direct CP violation. 



2 Search for K + ~ *-k + vv 
2.1 Experimental Design 

Since the expected branching fraction is « 10~ 10 , backgrounds are the major 
concern, consequently, the entire experiment and the analysis techniques are 
geared towards reducing backgrounds while maintaining a reasonable detection 
efficiency. 

The two major backgrounds are K + — > n + v (BR=64%) and K + — > 7r + 7r° 
(BR=21%). Other sources, e.g., K + — ► Tr°e + v e , if + — ► 7r°/i + ^ M , etc., are 
not as daunting. Fig. [I] shows the momentum spectrum of charged tracks in 
major K + decays. The experimental strategy is set by the characteristics of 
the backgrounds. 
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Figure 1 : Momentum spectrum of charged tracks produced in K + decays 

Since the two largest decay modes give off mono-energetic charged daugh- 
ters (Pjr = 205 MeV/c and = 236 MeV/c, respectively), we stop the inci- 
dent beam and let the K + decay at rest. If we maintain very good kinematic 
resolution, and search in a momentum region inside this range, we have one dis- 
criminant against these two background sources (the momentum spectrum in 
Fig. for the 7r + in K + — > 7r + vv decays has been enhanced by ~ 10 10 ). We also 
use the energy and range of the tracks as part of the kinematic discriminant. 

To further discriminate against muons, we detect the decay chain, 
7r + — » ji + —> e + . To do this, we require that the outgoing 7r + //i + stop in the 



detector and decay. The tt + will generate a /i + and an e + , whereas a muon 
will only give rise to an e + . 

Most K + decays produce a ir° in the final state, so we designed the ex- 
periment to have very good photon veto capabilities, e.g., the average rejection 
for the 7T° in K + — » 7r+7r° is about 10 6 . 

There are other sources, e.g. pions in the incident beam which scatter in 
the target and fall inside the signal region, charge exchange where a K + turns 
into a Kl, and the latter decays semi-leptonically. Such sources can be cut out 
by requiring that the outgoing pion be detected a few nano-seconds after the 
incident K + enters the target. This not only drastically reduces the possibility 
that an incoming pion scatters into the fiducial region, but also ensures that 
the incoming K + decays at rest. 

The E787 has been described in detail elsewhere © . Its main features 
are a degrader to slow the incident K +, s, so that they come to rest in an 
active target which is mainly composed of 413 5mm scintillating fibers. The 
target is surrounded by a drift chamber, and a range stack composed of 21 
layers of plastic scintillators. The outgoing daughter pions are required to stop 
in the range stack, where we detect the 7r — * fi — > e decay chain using 500 
MHz transient digitizers. Photons are detected with an extensive system of 
photon veto detectors, e.g., Pb-Scintillator calorimeter in the central region, 
Csl crystals in the endcap regions, etc. Track momentum is measured in the 
presence of a IT solenoidal magnetic field. 

2.2 "Online" improvements 

Since we first published the observation of this decay mode H) using data col- 
lected in 1995, we have taken more data during runs in 1996, 1997 and 1998. 
We have made many improvements; lowering the momentum of the K + beam, 
which results in a higher fraction of incident K + stopping in the target and a 
reduction in background hits from K + interactions in the degrader, improve- 
ments to the trigger, increased acceptance for detecting fi — ► e decays, etc. 

In Table [j], we present some details of the data runs in 1995-97. 

Even though the data runs were shorter in 1996 and 1997, we still collected 
more data as compared to the run in 1995. This was entirely due to the 
improvements outlined above. 



Table 1: 1995-91 dataset 





1995 


1996 


1997 


Total 


Length of run (weeks) 
Triggers (10 12 ) 
No. of evts. to tape 


25 
- 1.53 


17 
- 1.16 


8 

- 0.59 


~ 3.28 
~ 3.1 x 10 s 



3 Offline Analysis 

Since the backgrounds are orders of magnitude larger than the signal, the 
analysis has to be designed very carefully. Some of the main features are, 

(a) Blind Analysis The signal region is hidden (by inverting cuts) while 
cuts are developed and background levels are estimated. To avoid any bias in 
the background estimates, we develop cuts on 1/3 of the dataset and take the 
(actual) background level from the remaining 2 /3 of the data and extrapolate 
the latter to the full sensitivity. 

(b) "Bifurcated" Analysis In the analysis, we make an "a priori" identi- 
fication of the background sources and then develop at least two independent 
cuts for each source. In this manner, we can measure the rejection of each 
cut using the data itself. For instance, K + — > ir + ir° is rejected using kinematic 
cuts based on measuring the Range, Momentum and Energy of the outgoing tt + 
AND photon veto cuts which work on detecting the ir°. Similarly, K + — > n + v^ 
is rejected by applying kinematic cuts on the \x and by requiring the presence 
of the 7r + — > fi + —> e + decay chain. 

The goals for the analysis of the 1996-97 datasets and the re-analysis of the 
1995 dataset were (a) to increase rejection so as to keep the total background at 
the same level as in the original result D, (b) maintain or (possibly) increase 
detection efficiency, and (c) devise methods to increase background samples 
which would lead to a better understanding and a more precise estimation of 
the background. 

A lot of improvements were made in the analysis; tracking changes in 
the target, range stack and drift chamber improved the resolution in range 
and momentum, a better electron (from [i decay) finding algorithm was put in 
place, new cuts were devised to reject K + decays in flight, etc. For instance, 
the new analysis had a 30% larger rejection (with the same acceptance) against 
muon backgrounds. 



The background estimates for the 1995-97 dataset are presented in Table^J 



Table 2: Background Estimate for 1995-97 datasets 





Total 


K+ -> TT+TT 

K+ - 
l(&2)-bcam 
Charge Exchange 


0.020 ± 0.01 
0.030 ± 0.01 
0.020 ± 0.02 
0.010 ± 0.01 


Total 


0.08 ± 0.02 



The background estimate in the published result H) was 0.08 ±0.03, which 
implies that we increased rejection by a factor of about 2.2. At the same time, 
we were able to increase acceptance from 0.16% to 0.21% (an increase of almost 
30%). 
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Figure 2: Result for 1995-97 datasets 



In Fig.|[ we present the results from the 1995-97 datasets. The plot shows 
the data after all cuts, except those on Energy and Range of the 7r + , have been 



made. We have one event in the signal region. Using the total luminosity and 
detection efficiency, we obtain @) 



B(K+ -> tt+i/F) = l.5±H x 1(T 10 (2) 

Assuming unitarity of the CKM matrix and V c b, this result implies, 
0.002 < \V td \ < 0.04. 

4 Results on Radiative Decays 

4.1 A'+^7r+7r°7 

We have 19,836 events in this decay channel H) , which represents an increase 
of a factor of eight in statistics over the previous best measurement. 

The data is expressed in terms of a variable W, which is defined as, 

W 2 = (p • q)/n4+ x (p + • q)/m 2 + (3) 
= E 2 x (E^ - x cos 6^+ 7 )/(m 2 <+ x m^ + ) 

where, p, p_|- and q are 4-momenta for the K + , ir + and gamma, respectively. 
The structure dependent (DE) contribution is mainly at high values of W, as 
shown in Fig. |^. Low values of W are populated by Inner Bremsstrahlung (IB) 
decays. 
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Figure 3: (left) W spectrum of the signal events and best fits to IB+DE [solid 
curve] and IB alone [dashed curve]; (right) W spectrum normalized to the IB 
spectrum. 

We measure, 

BR(K+ -> 7r+7r°7; DE) = (4.72 ± 0.77) x 10~ 6 (55 < T..+ < 90MeV) (4) 



This result is roughly a factor of 4 lower than previous measurements. In 
addition, the interference, between IB and DE, is measured to be (— 0.4±1.6)%, 
and the ratio, DE/IB is measured to be (1.85 ± 0.30)%. 

The decay rate, corrected to full phase spacef], is now measured to be 
similar to that for K L : T(K+ -^tt+tt ^; DE) = (808 ± 132)s~ 1 vs. 
r(K° L ->TT + TT- r , DE) = (617 ± 18)s~ 1 . 

4.2 K+^/i+i/^ 

The structure dependent (SD) contribution has two components, due to the 
two polarizations of the outgoing photon. SD + peaks at higher E^ and E 7 and 
is easier to detect. We have roughly 2700 events in the region of interest, 
E 7 > 90MeV and E^ > l37MeV. 

In Figure [|, we present the data in terms of the angle between the /i + 
and the 7. It is clear from the fits that there is a large contribution from SD. 
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Fi gure 4: Distribution of cosO^y for E787 — ► fx^u^ candidates. A fit to 
IB alone is shown in (a), and a fit to both IB and SD + is shown in (b). 

The final result is based on a fit which includes contributions from SD~, 
IB and interference terms, 

B{SD+) = (1.33 ±0.12 ±0.18) x 10~ 5 (5) 

The sum of form factors is measured to be, \Fy + Fa\ = 0.165 ± 0.007 ± 
0.011. This can be compared to the 0(p 4 ) xPT calculation 0) of Fy + Fa = 
0.137 ± 0.006 and B(SD+)=9.22x 10~ 6 . We also measure -0.04 < F v - F A < 
0.24 at 90% CL; this is to compared with an expectation 0) of (0.052 ± 0.006). 

^^This correction assumes that the form factor has no energy dependence. 



5 Future Plans 



5.1 1998 dataset 



As mentioned in Section 2.2, we also took data in 1998. The sensitivity of this 
data set is expected to be approximately the same as that of the combined 
1995-97 datasets. This result is expected in 2001. 



5.2 Phase Space below K^2 peak 

In the analysis presented here, we have searched for the signal in the region 
above the and below the momentum peaks. However, because of the 
V-A nature of the the decay, there is also a large phase space below the K W 2 
peak, as shown in Fig.|l|. Some of the advantages of looking in this region are, 
phase space is larger than above K^2, less tt + absorption losses, etc., leading 
to a possibility of x3 increase in acceptance. 

However, the background due to K n 2 is very large, since the 7r + can elas- 
tically scatter in a fiber hit by the incoming K + , lose energy and fall inside 
the signal region. Since this is a real pion, we cannot use the kinematic dis- 
criminant to reject this background. This scenario occurs when the daughter 
pions travel along the beam axis. Since photon veto capability is weak in this 
angular region, it is more likely to miss the two decay photons. One has to 
look for evidence of the pion scatter in the kaon fiber, as a way to discriminate 
against this background. 

We made good progress in studying this region using the 1996 dataset. 
Searching in the momentum region, 140 MeV/c < P^ < 190 MeV/c, prelimi- 
nary studies showed that Ssm/B was ~ 1:5, where Ssm is the expected Standard 
Model signal. We expect this ratio to improve to ~ 1:3 when we extrapolate 
to 1995-98 datasets, and to ~ 1:1 in experiment E949. 



5.3 Experiment E949 

Experiment E949 is an upgrade of E787 and is scheduled to run in 2002-2003, 
with an engineering run scheduled for 2001. 

Many subsystems in E787 have been upgraded, e.g., Photon veto has been 
improved by adding more detectors in the central region, converting the current 
degrader into an active degrader, adding a photon veto detector downstream 
of the target, etc., the Range Stack has been improved by replacing some of the 



scintillator layers, Trigger and DAQ systems have been upgraded, etc. These 
improvements are expected to increase sensitivity by about a factor of 3 (over 
the result based on the 1995 dataset H)). 

In addition, the accelerator is scheduled to deliver more protons per spill, 
increase the duty cycle and reduce the K + momentum, leading to an improve- 
ment of a factor of 2.2 in sensitivity. 

We plan to reoptimize our analysis for these higher rates and also hope 
to gain sensitivity from looking in the phase space below the K^2 peak, as 



described in Section 5.2. Each of these two sources is expected to yield an 
increase in sensitivity of a factor of 2. However, we take a conservative approach 
and take only one factor of 2 in our projections. 

Factoring in the increase in running time (~ 6000 hours), we estimate 
that our sensitivity will be a factor of 50 higher than the result based on 
the 1995 dataset u) , which implies a Single Event Sensitivity (SES) of about 
8 x 10~ 12 . This means that we expect to observe 10 Standard Model events 
in E949, allowing us to measure V t d with a precision of about 23%, which is 
extremely competitive with the current precision of about 21% obtained from 
B d B d mixing. 



6 Conclusions 



Using data collected in 1995-97, we have measured , 

B(K+ ir+W) = 1.5+H x 10- 10 (6) 

We have also made high statistics measurements of the structure depen- 
dent parts of i^ + ^Tr + 7r°7 and K + ^n + v^ . 

Using all the data collected with the E787 detector, we expect to improve 
the SES for K+^tt+W to - 0.7 x lO" 10 . 

Using the data which we expect to collect in experiment E949, we hope 
to observe O(10) K + -^>tt + vT> Standard Model events. This will allow us to 
measure V^with a precision of about 23%. 

In Fig.^, we present the improvement in SES. We also show the expecta- 
tions from the 1995-98 datasets (denoted as E787 final), from experiment E949, 
and from a proposed experiment (CKM) at FermilabB). 
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Figure 5: Single Event Sensitivity 



References 



1. A. Buras and R. Fleischer, Los Alamos e-archive, hep-ph/970437(: 



2. T.K. Komatsubara et. at, Nucl. Instrum. Methods Phys. Res., Sect A 404, 
315 (1198) and references within 

3. S. Adler et. ai, Phys. Rev. Lett. 79, 2204 (1997) 

4. S. Adler et. ai, Phys. Rev. Lett. 84, 3768 (2000) 

5. S. Adler et.al Phys. Rev. Lett. 85, 4856 (2000) 

6. S. Adler et. al. Phys. Rev. Lett. 85, 2256 (2000) 

7. J. Bijnens et al, Nucl.Phys. B396, 81 (1993) 

8. See talk by Hogan Nguyen at this conference 



